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Abstract. Measurement of the transverse spin dependence of the di-jet opening angle in pp col- 
lisions at s/s = 200 GeV has been performed by the STAR collaboration. An analyzing power 
consistent with zero has been observed over a broad range in pseudorapidity sum of the two jets 
with respect to the polarized beam direction. A non-zero (Sivers) correlation between transverse 
momentum direction of partons in the initial state and transverse spin orientation of the parent pro- 
ton has been previously observed in semi-inclusive deep inelastic scattering (SIDIS). The present 
measurements are much smaller than deduced from predictions made for STAR di-jets based on 
non-zero quark Sivers functions deduced from SIDIS, and furthermore indicate that gluon Sivers 
asymmetries are comparably small. 
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Experimental determination of the contribution from orbital angular momentum of 
the proton's constituents to its spin is difficult. A possible manifestation of orbital 
momentum are Sivers correlations [Q]] between the transverse spin (St) direction of a 

polarized proton and transverse momentum (k^) directional preferences of unpolarized 
partons carrying longitudinal momentum fraction x of the proton: 

, P-fk^xSr) 

J- (1) 

2 \P\-\k L \-\S T \ 

where P is the momentum vector of the polarized proton and A N f(x,k ± ) is called the 
Sivers function. 

Recently measurement of a non-zero Sivers function for tz + in SIDIS has been 
reported by the HERMES collaboration [2]. The same group found an asymmetry 
consistent with zero for n~ , leading to an interpretation [II that the Sivers functions are 
opposite in sign and different in magnitude for u vs. d quarks. Boer and Vogelsang [0] 
suggested that the Sivers effect would be manifested as a spin-dependent azimuthal 
side preference for almost back-to-back di-jets produced in proton collisions involving 
transversely polarized beam. 

In this paper we report Sivers single spin asymmetry (An) measurements for di-jet 
production in kinematic regions that allow us to distinguish quark from gluon Sivers 
functions. We compare our results to predictions [3] derived from SIDIS, in order to 
judge the universality of the extracted quark Sivers functions. 




FIGURE 1. Geometry of the STAR detector. 



The di-jet measurement has been performed in the 2006 run using the STAR detec- 
tor at RHIC JHt]. A 3 week- long run resulted in a sampled luminosity of 1.1 pb _1 . Two 
polarized beams collided in the center of the STAR detector, as shown in Fig. [TJ Elec- 
tromagnetic calorimeters (EMC) covering pseudorapidity range r\ E [— 1 , +2] and full 
azimuthal angle have been used to trigger on di-jets. A dedicated software trigger 
accepted events with two EMC patches covering 0.6 x 0.6 in At] x A0, separated az- 
imuthally by at least 60°, and containing at least 3.3 and 3.6 GeV of transverse energy 
(Et), respectively. Additionally, the hardware trigger required Ej > 14 GeV summed 
over the whole EMC and a coincidence between 2 forward segmented scintillators 
(BBC) placed at 3.3 < |t] \ < 5.0 on either side of the interaction point. In this analy- 
sis we have used the on-line values of T] and for both jets obtained from iSj-weighted 
centroids of the EMC tower locations included in the jet patches. The results reported 
here include only those events where both detected jets had EMC Ej < 8.0 GeV within 
their respective 0.6 x 0.6 At] x A0 patch. 

The azimuthal correlation between the two jets is shown in Fig. [2k, where back- 
to-back di-jets dominate the spectrum. For the Sivers analysis, we define the "signed" 
azimuthal di-jet opening angle £ (Fig. [2b), chosen to be > 180° when cos(j)b > and 
< 180° otherwise, where § b is the di-jet bisector angle measured with respect to STAR'S 
+x-axis (see Fig. [I]). The spectrum of average r,l ^ ri2 is shown in Fig. [2h. The most 
forward measured di-jets when T]i + T]2 > 2 involve preferentially quark (gluon) from 
the polarized beam heading toward +z (-z) direction. 

The vertical beam spin orientation alternated for each successive bunch of one beam 
and each pair of bunches of the other. An integrated Sivers was extracted for each 
beam from spin-sorted ^-spectra using a cross-ratio in which "left" and "right" yields 
were deduced from £ > 180° and £ < 180°, respectively, for the +z beam and the op- 
posite assignment was used for the — z beam. On-line values of beam polarization have 
been used and a systematic normalization error of ±20% is assigned to until final 
beam polarization values are known. The cross-ratio method does not require indepen- 
dent luminosity monitors if the detector has approximately equal left-right efficiency, as 
is true for this measurement. 

The jet direction deduced from its electromagnetic component only can deviate from 
the direction of the parent parton due to fragmentation processes and finite detector 
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FIGURE 2. Characteristic distribution of events that pass the STAR di-jet trigger: a) azimuthal correla- 
tion, b) signed opening angle and c) averaged pseudorapidity, all from trigger-level EMC information 
only. 



resolution. Simulations with the Pythia 6.205 event generator and full GEANT modeling 
of the STAR detector have been used to estimate this azimuthal smearing to be o ~ 6.3°. 
This ^-smearing is small compared to the observed overall £ width of ~ 20°, driven 
by the intrinsic k smearing of the partons inside the proton. However, this 0-smearing 
does cause a dilution of the measured An by 10% of the input parton-level asymmetry, 
as determined from a toy model simulation. 

The toy model assumed a binary parton collision with an initial Gaussian k 1 - dis- 
tribution with added exponential tails. To simulate the Sivers effect it assumed a small 

offset in the centroid of the model k L distribution, in the direction St x P, using the 
randomly assigned spin direction for each beam. Azimuthal directions of the outgoing 
partons were smeared with a Gaussian of o = 6.3° to mimic fragmentation and detec- 
tion processes. The toy model parameters were tuned to reproduce the observed jet Ej- 
(from full jet reconstruction including charge-particle tracks) and ^-spectra. 

Fig. [3] shows the measured analyzing power for both beams (heading toward ±z) 
compared to predictions [3] as function of r\\ + r\2 integrated over the full STAR EMC 
r\ acceptance and for |£ — 1 80° | < 68°. The full analyzed sample yields 2.6 x 10 6 di- 
jets with average A^ z values consistent with zero for both beams, within statistical 
uncertainties of ±0.002. Based on the toy model these results probe a Sivers (H) 
preference as small as ±3 MeV, or ±0.2 % of the inferred rms width of the k X j 
distributions. 

The systematic error bands in Fig. [3] arise from two sources combined in quadrature. 
An error of ±20% of the statistical uncertainty on each bin allows for multi-jet event 
contribution to the £ distribution wings (see Fig. |2b). In addition, we include a model- 
dependent uncertainty on the dilution factor (0.90 ±0.05) deduced from the toy model. 
The 20% uncertainty on the on-line beam polarization values is not included in the error 
band. 

The predictions fl] in Fig. El are based on two models (labeled VY1 and VY2) of 
the u and d quark Sivers functions fitted to HERMES SIDIS results [2]. The models 
differ in tying the d-quark Sivers function shape to either the u{x) (VY1) or d{x) 
(VY2) unpolarized parton distribution function. Both calculations assume zero gluon 
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FIGURE 3. Compai-ison of STAR results for A^ z (left) and A^ z (right) as function of 77 [ + 772 with 
predictions H], based on two models of quark Sivers functions deduced from fits to HERMES SIDS 
results Jl. 



Sivers function. The predictions are integrated over a jet pj range (5-10 GeV/c) well 
matched to our measurements, and further over the STAR detector acceptance within 
each r/i + 772 bin. In these calculations, di-jet production was treated like Drell-Yan di- 
lepton production, with gauge link factors for initial-state interactions only [3]. We have 
reversed the sign of the predictions to adhere to the Madison convention used for our 
measured asymmetries. 

In the region 771+ 772 > 2, A^ z (AZ Z ) preferentially probes quark (gluon) Sivers 
functions. The predicted A^ z values in this region reflect the sizable HERMES Sivers 
asymmetries, while A^ z ~ because gluon Sivers effects have been ignored. 

All measured di-jet asymmetries are within two standard deviations of zero, suggest- 
ing that both quark and gluon Sivers effects are considerably smaller for di-jet produc- 
tion in pp collisions than are the quark Sivers effects inferred from the HERMES results. 

The experiment-theory discrepancies in Fig. [3] might reflect a breakdown of pQCD 
factorization for the soft values probed for back-to-back di-jets (near £ = 180°). It 
remains an open challenge for theory to account simultaneously for the present di-jet 
results and for the sizable transverse single-spin asymmetries seen in SIDIS [2] and in 
forward pp —> 7t°X ((J. 
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